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The CALICE collaboration develops prototypes for highly granular calorimeters for 
detectors at a future linear electron positron collider. The highly granular electromag- 
netic calorimeter prototype was tested in particle beams. We present the study of the 
interactions of hadrons in this prototype. 
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i— i 1 Introduction 
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The next machine after the LHC will be a linear electron positron collider at the TeV scale. 

A This machine will allow high precision measurements to extend the scientific results of the 

LHC. Currently the most advanced proposal is the International Linear Collider (ILC) with 

*—J a centre-of-mass energy of up to 500 GeV. An alternative at higher center-of-mass energies 

than the ILC is the Compact Linear Collider (CLIC). 
• i-H The e + /e~ collisions will lead to multi-jet final states. Their reconstruction will be 

based on so-called particle flow algorithms (PFA) [TJ|2]. The goal is to reconstruct every 
single particle of the final state which requires the combination of information from each 
sub-detector. Particularly a perfect association of the signals in the tracking systems with 
those in the calorimeters is required. Charged particles are measured with the tracking 
system, while neutral particles can only be measured in the calorimeters. The particle flow 

^ technique requires a specific detector design. A limiting factor of the particle flow, called 

confusion, is due to overlapping showers. To separate the particles, calorimeters need to be 
of high granularity. They have to be very compact since they will be embedded within the 
magnetic coil of the detectors. 

The CALICE collaboration [3] designs and studies electromagnetic and hadronic calorime- 
ters for experiments at the ILC. All detectors are developed with a common approach to 
obtain ultra granular calorimeters optimised to particle flow algorithm. The prototypes 
are tested in combined beam tests in order to study several technologies and the associa- 
tion of electromagnetic and hadronic calorimeters. These beam tests allow for developing 
. ^h reconstruction algorithms and for the validation of simulations. 

>< 

2 Silicon tungsten electromagnetic calorimeter 

The silicon tungsten electromagnetic calorimeter (ECAL) is built following particle flow 
algorithm requirements and based on a sampling design. The physics prototype of the 
ECAL, used in this study, is composed of 30 layers of silicon as active material, alternated 
with tungsten as absorber material. Tungsten has a short radiation length and a small 
Moliere radius which gives compact showers and allows for the efficient separation of close 
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particles. The silicon allows for a thin and easily segmented readout detection system suited 
for high granularity. The silicon wafers are segmented in 1 x 1 cm 2 pads. Each wafer consists 
of a square of 6 x 6 pixels and each layer is a matrix of 3 x 3 wafers resulting in an active 
zone of 18 x 18 cm 2 giving a total of around 10000 readout channels. As shown on Figure [I] 
the printed circuit boards (PCB) are mounted two by two in an elementary detection unit 
called slab. The silicon wafers are located on each side of a H-shaped tungsten supporting 
structure. These slabs are shielded by aluminum foils and inserted in the a carbon fiber 
composite mechanical structure. Three different longitudinal samplings were used, finer in 
the first ten layers (1.4 mm plates of tungsten), 2.8 mm plates in the ten intermediate layers, 
and 4.2 mm in the last ten layers 4 . 
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Figure 1: Pictures of the ECAL prototype. The top figure shows the structure of a single 
detector slab while the bottom figure shows an overview of the whole prototype. 

During the years 2005-2011, the calorimeter was tested with various beams: electrons, 
positrons, muons, pions and protons in the momentum range between 1 and 180 GeV/c. 
Beam campaigns were conducted at DESY, CERN and FNAL. The main goal of these beam 
tests was to demonstrate the principle of highly granular calorimeters. Using data collected 
with electron beams at energies from 6 GeV to 45 GeV, the linearity and resolution of this 
calorimeter to electrons were estimated |5 . The linearity is better than 1% while the relative 
energy resolution is estimated to be cr E /E = (16.53 ± 0.14(stat) ± 0.4(syst))/ v / ^(GeV) © 
(1.07±0.07(stat) ±0.1(syst)) (%). To compare the test beam data with the simulation, the 
detector has been simulated within the Geant4-based MOKK A simulation framework 6 7 . 

In addition to standard beam test analysis (calibration, noise and stability), the high 
granularity allow to study hadronic interactions in the detector and to validate hadronic 
model implementations in Geant4 simulations. As hadronic interactions are poorly under- 
stood, the high granularity of the calorimeter offers unprecedented information to study 
hadronic interactions. Several model combinations called physics lists are proposed in 
Geant4. Our studies allow for a detailed testing of these simulations. A study of inter- 
actions of pions in the physics prototype has been published in M for energies between 
8 GeV and 80 GeV. The work presented in this article concentrates rather on energy range 
2 GeV -10 GeV. 
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3 Results 

The high granularity permits detailed view into hadronic showers. In order to exploit the 
imaging of hadronic showers, we have to develop new analysis methods like particle tracking 
or interaction vertex localisation in the calorimeter. 

The high granularity of the calorimeter allows the tracking of particles as they pass 
through the detector and the use of imaging processing techniques. For instance, the Hough 
transform technique has been tested to find tracks in the calorimeter |9 . This technique has 
been applied to reconstruct a muon track near to a 30 GeV electromagnetic shower. The 
Figure [2] shows the reconstruction efficiency as a function of the distance between the track 
and the shower axis. The efficiency reaches 100% for distance of 25 mm. This separation 
efficiency is essential for particle flow algorithm, especially to separate charged and neutral 
particles which are not detected by trackers. 
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Figure 2: Efficiency of MIP detection as a function of distance from electron shower 



Figure [3] illustrates the potential of the high granularity. We distinguish two different 
and recognizable topologies. Events can be classified as a function of their topologies into 
4 classes. MIP and elastic scattering are considered as non interactive events. The two 



distributions presented in Figure [3 
(top) and "Fireball" (bottom) MO 



jshow the two types of interacting events called "Point like" 

.m- 

In Figure [4j we compare the relative frequency of each event class between data and 
some Geant4 physics lists. The frequency of "Fireball" events in the data is about 55% over 
the whole energy range. After delta-rays correction, the frequency of "Pointlike" events is 
about 4% at small energies and tend towards zero at higher energies. The frequency of both 
classes for inelastic events is well reproduced by all physics lists which confirms that the 
total inelastic cross sections are well implemented into the physics lists. Cross sections of 
non-interacting events are also well modeled by Geant4 [To] 
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3.1 Transverse profiles 

The shower radius is defined by: 



Rf 
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Figure 3: Examples of hadronic interactions in the ECAL. "FireBalF event on the top and 
"PointLike" event on the bottom. 
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and the same for y. The transversal profile is an important observable affects the overlap of 
shower thus the confusion. 

For the calculation of the observables, only hits in the interaction layer and all subsequent 
layers are taken into account. In order to define in the same way a measure of the radius 
of non interacting events, i.e. interactions where no interaction point could be found, the 



width Re is calculated by summing over all hits in the ECAL 10 11 



The left part of Figure [5] shows the comparison of the transverse profile between data 
and simulation at energies from 2 GeV up to 10 GeV. In this plot we use the QGSP_BERT 
physic list. We distinguish 2 maxima at 5 mm and 20 mm. Up to 6 GeV, data are well 
reproduced by the simulation except for radius larger than 35 mm. At 8 GeV the simulation 
overestimate the data in the transition region between the two peaks. The situation is 
comparable at 10 GeV 

In the right part of Figure |5j the contribution of each event type is separated for 8 GeV 
events. We see that the first maximum in Figure [5] left is created by MIP events and the 
second one by "Fireball" events. The maxima of "Pointlike" and "Scattered" events are 



located in the transition region. The four event types show the same disagreement 10 11 



3.2 Longitudinal profiles 

The longitudinal profile is given as a function of pseudolayers in order to account for the 
different sampling fractions in the ECAL. In the first module, pseudolayers are equivalent 
to real layers. In the second module, each layer is subdivided in two pseudolayers and 
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Figure 4: Frequency of event types and comparison with several physics lists as function of 
energies from 2 GeV to 10 GeV. The top figure shows the rates of interactive classes while 
bottom figure shows the rates of non-interactive classes. 



in the third module, layers are subdivided into three pseudolayers. The energy is linearly 
interpolated between the pseudolayers. In case of interacting events, the longitudinal shower 
profile starts from the interaction layer. In case of non-interacting events, the longitudinal 
profile is calculated from the first detector layer 
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In Figure |6j the test beam data are shown in red and the simulation in black. For 
simulated events, we separated the contributions of the secondary particles. The blue lines 
are the contributions from electrons and positrons. The green lines are the contributions 
from protons. The red lines are the contributions of pions. Other particles are in purple. For 
"Fireball" events, none of the models give a satisfactory description of the data. For instance, 
the top distribution shows significant disagreement between data and the QGSP_BIC model. 
In the bottom part of Figure [6j we present longitudinal profile of "Pointlike" events and the 
comparison with QGSP-BERT model. This class of events is better reproduce by the Bertini 
model 
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4 Conclusion and outlook 

This study investigates in depth the interactions of pions in an energy range between 2 
and 10 GeV. It demonstrates the outstanding potential of the CALICE silicon tungsten 
electromagnetic calorimeter to resolve details of hadronic cascades. The data obtained in 
test beams are compared with different physics lists provided by Geant4. Visible differences 
between models appear at this level of detail. The study on shower decomposition is the first 
step toward the use of sophisticated algorithms which may be developed in collaboration 
with experts from applied mathematics. To succeed the physics prototype, a technological 
prototype is in development which will feature a four times higher granularity. 12 
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Figure 5: Left: Lateral distributions of the hadronic shower energy at energies from 2 GeV 
to 10 GeV. For each distribution: the top view features the comparison between test beam 
data in red and the simulation with QGSP-BERT in black. The bottom view shows the 
ratio of test beam data over simulation. Right: Lateral distributions of the hadronic shower 
for each event classes at 8 GeV. For each distribution: the top view features the comparison 
between test beam data in red and the simulation with QGSP-BERT in black. The bottom 
view shows the ratio of test beam data over simulation. 
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Figure 6: Total longitudinal profiles for the two interacting event types: comparison between 
test beam data in red and simulation in black at 2 GeV. For simulated events, we separated 
the contributions of the secondary particles. The blue lines are the contributions from 
electrons and positrons. The green lines are the contributions from protons. The red lines 
are the contributions of pions. Other particles are in violet. The top figure shows the 
comparison between test beam data and the QGSP-BIC physic list for "Fireball" events 
while the bottom figure shows the comparison between test beam data and the QGSPJBERT 
physic list for "Pointlike" events. 
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